Pion and kaon production in central Pb+Pb collisions at 20A and 30^4 GeV: 
Evidence for the onset of deconfinement 
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Results on charged pion and kaon production in central Pb+Pb collisions at 2QA and 30^ GeV 
are presented and compared to data at lower and higher energies. A rapid change of the energy 
dependence is observed around 30j4 GeV for the yields of pions and kaons as well as for the shape 
of the transverse mass spectra. The change is compatible with the prediction that the threshold for 
production of a state of deconfined matter at the early stage of the collisions is located at low SPS 
energies. 

PACS numbers: 25.75.-q 



INTRODUCTION 



The advent of the quark model of hadrons and the development of quantum chromodynamics naturally led to the 
question whether strongly interacting matter exists in different phases and, if so, of which nature the transitions 
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between these phases are. In particular, it is commonly believed that a gas of hadrons will undergo a transition to 
a state of quasi-free quarks and gluons, the Quark Gluon Plasma (QGP) when its temperature exceeds a critical 
value These questions have motivated a broad experimental program of nucleus- nucleus collisions to study the 
properties of strongly interacting matter at extreme densities and temperatures. 

Several signatures of the formation of a transient QGP state during the early stage of a nucleus-nucleus collision 
at high energies were proposed in the past However, the validity of these signatures has come under renewed 
scrutiny. For this reason the NA49 Collaboration at the CERN SPS has searched over the past years for signs of the 
onset of QGP creation in the energy dependence of hadron production. This search was motivated by the prediction 
of [1, m 0] that the onset of deconfinement should lead to a steepening of the increase of the pion yield with collision 
energy and to a sharp maximum in the energy dependence of the strangeness to pion ratio. The onset was expected 
to occur at approximately SO A GcV @. 

The NA49 energy scan program at the CERN SPS started with two runs where data on central Pb-|-Pb collisions 
at 40A and 80A GeV were recorded in 1999 and 2000. Data at the top SPS energy of 158A GeV had already been 
taken in previous SPS runs. The analysis of these runs was published in [7] and the results confirmed the predictions 
of 0]. This finding motivated an extension of the energy scan to the lower SPS energies of 30 A and 20^ GeV which 
was completed in 2002. In this letter we report final results on charged kaon and pion production in central Pb-|-Pb 
collisions at 20A and 30^ GeV beam energy. The present measurements are combined with those from Q and 
compared to available model calculations, together with lower and higher energy data from AGS and RHIC. 



DETECTOR AND DATA ANALYSIS 

The NA49 experimental set-up Q consists of four large volume Time Projection Chambers (TPC). Two of these 
(VTPC) are placed in the field of two super-conducting dipole magnets. The other two TPCs (MTPC) are positioned 
downstream of the magnets and are optimized for high precision measurements of the ionization energy loss dE/dx 
with a resolution of about 4%. This dE/dx measurement provides particle identification which is complemented by 
a measurement of the time-of-flight (TOE) with a resolution of about 60 ps in two TOE detector arrays positioned 
downstream of the MTPCs. At each incident energy the TOE acceptance for kaons was kept at mid-rapidity by 
lowering in proportion to the beam energy the nominal 158^4 GeV field settings of about 1.5 (upstream magnet) and 
1.1 T (downstream magnet). In both the 30 A and 20 A GeV runs a thin lead foil target of 224 mg/cm^, corresponding 
to about 1% of interaction length for Pb ions, was positioned 80 cm upstream from the first TPC. A trigger based on 
a measurement of the energy deposit of projectile spectator nucleons in a downstream calorimeter selected the most 
central 7.2% of the Pb-|-Pb collisions. The corresponding mean number of wounded nucleons A^w 0] is calculated 
using the Eritiof model [l3| to be = 349 ± l(stat) ± 5(sys). At each energy about 3.5 x 10^ events were recorded. 

Results on kaon and pion production were obtained in a multi-step analysis procedure which involves charged 
track reconstruction, particle identification and corrections to account for background contributions as well as for 
acceptance and efficiency losses. To reduce the systematic errors, the analysis has been restricted to regions of 
phase space where background and efficiency corrections are small and approximately uniform. To minimize tracking 
efficiency corrections, only tracks within an azimuthal angle wedge of ±30° with respect to the horizontal plane were 
used. 

The spectra of tt^ and at mid-rapidity are obtained using the combined dE/dx and TOE information as 
described in [lH (TOE + dE/dx analysis). The pion spectra were corrected for the contribution from the weak 
decays of the strange particles A, and Kg as well as for ^ contamination (charged products of weak decays 

can be reconstructed as coming from the event interaction vertex, while muons cannot be separated from pions in 
NA49 with dE/dx and TOE measurements for momenta above about 1 GeV/c). The correction factors were obtained 
from a GEANT simulation of the NA49 detector using as input the hadron distributions from the VENUS model [13] ■ 
The model distributions were tuned to reproduce the measured yields of A [l^l and Kg where the latter were taken 
to be the average of the yields presented in this paper. The total background correction to the tt" (tt"^) yields is 
found to be about 8 (6)%. 

Raw and K~ yields at forward rapidities were extracted from fits of the dE/dx distributions in narrow bins of 
momentum and transverse momentum. The fitted function parametrizes contributions from e+ , n'^ , K"*" , protons and 
deuterons to the dE/dx spectra of positively charged particles and corresponding contributions to those of negatively 
charged particles. Eor pions the dE/dx method does not provide sufficiently accurate identification in the low pT 
region near mid-rapidity. Therefore to obtain the raw 7r~ spectra in the full forward hemisphere of the reactions 
yields of all negatively charged particles were determined as a function of rapidity (calculated assuming the pion 
mass) and transverse momentum px- The contamination by K", p and e^ from the interaction vertex as well as 
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non-vertex hadrons originating from strange particle decays and secondary interactions was subtracted. The total 
correction amounts to 20-25% and was calculated using the VENUS/GEANT simulation. The above procedure can 
not be applied for positively charged hadrons due to the large contribution of kaons and protons. 

The resulting background subtracted spectra were corrected for geometrical acceptance, losses due to in-flight 
decays, inefhciencies of the tracking algorithms (about 5%) and quality cuts. More details on the correction procedure 
can be found in Q. 

Systematic errors were estimated by comparing results obtained with different detectors (TPC, TOF) and by varying 
cuts and correction strategies. Uncertainties in the parameters of the AE /Ax fits lead to asymmetric systematic errors 
on the kaon yields. The systematic errors were estimated to be 5-10% and are explicitly given in Table II. Note that 
to a large extent the systematic uncertainties presented here are common to those reported in [3] since there the same 
experimental procedure was used. 



RESULTS AT 20A AND 30A GEV 



In Fig. [T] are shown the transverse mass spectra of tt^ mesons near mid-rapidity < y < 0.2 together with those 




FIG. 1: Transverse mass spectra of tt" and tt"*" mesons at < j/ < 0.2 (left) and those of K"*" and K~ mesons at \y\ < 0.1 
(right) produced in central Pb-|-Pb collisions at 20^4 and 30A GeV. The lines are fits of Eq. ^ to the spectra in the interval 
0.2 < — m < 0.7 GeV. The statistical errors are smaller than the symbol size. The systematic errors are 5% in the region 
used for the fit and reach 10% at the low and high ends of the mx spectra. 



of at |y| < 0.1 obtained from the TOF + AE/Ax analysis of 2QA and 30A GeV central Pb4-Pb collisions. Here 
the transverse mass is defined by = p\ -'r ™^ with m the rest mass of the particle, and y denotes the rapidity of 
a particle in the collision center-of-mass system. The full lines in Fig. [T] indicate a fit of the function 

d^'^ = C cxp f (1) 
mTdrnxdy V T / 

to the data in the range 0.2 < tot — m < 0.7 GeV. The values obtained for the inverse slope parameter T are given 
in Table U 

The kaon spectra are well described by the fit, while the low and high mx regions of the pion spectra outside of 
the fitted region are above the extrapolated fitted line. 

The rapidity distributions An /Ay of vr^, and mesons are plotted in Fig.[2l These distributions were obtained 
by summing the measured mx spectra and using the fitted exponential function Eq. ([1]) to extrapolate to full m^. 
For kaons the corresponding corrections are below 10% for most of the rapidity bins except for the first two and the 
last rapidity bin from the AE/Ax only analysis in which not more than about 50% of the yield is in the measured 
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FIG. 2: Rapidity distributions of it~ (left), K~ (middle) and K"*" mesons (right) produced in central Pb+Pb collisions at 20^4 
(lower curves) and 30A GeV(upper curves). Squares and circles show the results of an analysis based on d_E/da; only whereas 
triangles are obtained from a TOF + di5/da; analysis. The closed symbols indicate measured points while the open points are 
reflected with respect to mid-rapidity. The lines indicate flts of Eq. ([2)) to the data. The errors, which are often smaller than 
the symbol size, are statistical only. The systematic errors on the measurements are about ±5%. 



region. The corrections for pions are negligible. The rapidity spectra were parameterized by the sum of two Gaussian 
distributions positioned symmetrically with respect to mid-rapidity 0] 
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where (n), a and yo are fit parameters. The results of the fits are indicated by the full lines in Fig. [2] and the values 
of the parameters are given in Table [D Also given in this Table are the mean multiplicities which were obtained by 
integration of the fitted curves and the mid-rapidity yields that were taken to be the maxima of the curves. The 



TABLE I: In the columns on the left are listed the inverse slope parameters T obtained from a fit of Eq. ([1} to the m-x spectra 
at mid-rapidity, together with the parameters a and j/o from a fit of Eq. Q to the vr" and rapidity spectra. Only statistical 
errors are given. The values of a and i/o are significantly correlated. In the columns on the right are listed the full phase space 
yields (vr^) and (K*) together with dn/dy for vr^ and production at mid-rapidity. The first error is statistical, the second 
systematic. Note that (tt^) is not directly measured (see text). 



Parameter 


20A GeV 


30A GeV 


Yield 


20A GeV 


30A GeV 


T{tv+) (MeV) 


167 ±2 


175 ±2 


(tt-) 


221 ± 1 ± 11 


274 ± 1 ± 14 


T(7r-) (MeV) 


160 ±2 


169 ±2 


(7r+) 


190 ± 1 ± 9 


241 ± 1 ± 12 


r(K+) (MeV) 


219 ±5 


232 ± 5 


(K+) 


40.7 ±0.7 ±2.2 


52.9 ± 0.9tli 


r(K~) (MeV) 


193 ±9 


230 ±7 


(K-) 


10.3 ±0.1 ±0.2 


16.0 ±0.2 ±0.4 




0.837 ± 0.007 


0.885 ± 0.007 


dn/dy(7r^) 


84.8 ±0.4 ±4.2 


96.5 ±0.5 ±4.8 


a{K+) 


0.601 ±0.012 


0.722 ±0.026 


dn/dy(7r"'') 


72.9 ±0.3 ±3.6 


83.0 ±0.4 ±4.2 


a(K-) 


0.642 ± 0.035 


0.710 ± 0.032 


dn/dy{K-) 


5.58 ±0.07 ±0.11 


7.8 ±0.1 ±0.2 


2/0 (^") 


0.557 ± 0.009 


0.624 ±0.009 


dn/dy(K+) 


16.4 ±0.6 ±0.4 


21.2 ± O.Sto.l 


2/o(K+) 


0.606 ±0.014 


0.578 ± 0.030 








J/o(K-) 


0.34 ±0.06 


0.37 ±0.05 









mean multiplicity of 7r+ mesons and their mid-rapidity yields were calculated by scaling (tt^) with the tt^/tt" ratio 
measured at mid-rapidity in the TOF + dE/dx analysis. These ratios were found to be (0.86,0.88,0.90,0.91,0.93) 
for the data measured at (20, 30, 40, 80, 158)A GeV. The ratios at the latter three energies were used to recalculate 
the published values of the (vr"'') and the tt"*" mid-rapidity yields at these energies. The recalculated 7r+ multiplicities 
differ by several percent from the ones published in 0] . 
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REVIEW OF THE ENERGY DEPENDENCE 



In this section, the new results on tt and K production at 20^1 and 30^ GeV will be discussed together with 
published measurements at lower (AGS) and higher (SPS, RHIC) energies and compared to the corresponding data 
from p + p(p) interactions. Model calculations, which are shown by the curves in the figures below, will be discussed 
in the next section. 

In Fig. [3] is shown the mean pion multiplicity (tt) = 1.5 ((tt^) + (tt^)) per wounded nucleon {N„) as a function of 




FIG. 3: Left: Energy dependence of the mean pion multiplicity per wounded nucleon measured in central Pb+Pb and 
Au+Au [Tsl. [T^ collisions (full symbols), compared to the corresponding results from p + p(p) reactions (open circles). Right: 
Energy dependence of the difference between the measured mean pion multiplicity per wounded nucleon and a parametrization 
(see text) of the p + p data. The meaning of the full and open symbols is the same as in the left-hand plot. The lines show 
various model predictions described in the text. 

the collision energy, expressed by Fermi's measure [l^ 

n 1/4 



F = 



(y/SNN - 2mN) 



where ^/snn is the centre-of-mass energy per nucleon-nucleon pair and rrijv the nucleon rest mass. The use of F in 
Fig. [3] is motivated by model expectations, as will be discussed below. 

In Fig. [3] the measurements by NA49 are compared to results from the AGS ^iB'l and RHIC fi^ on central nucleus- 
nucleus collisions. The lowest point at the top RHIC energy is obtained from the BRAHMS measurements of identified 
hadron spectra. The middle and highest points are results from, respectively, BRAHMS and PHOBOS on the charged 
hadron mean multiplicity corrected for non-pion contributions by use of the BRAHMS identified hadron yields. The 
RHIC points at the lower energies are from PHOBOS, again corrected for the non-pion contribution. The results 
from p + p(p) interactions are shown by the open symbols. For ease of comparison these data were fitted to the 
parameterization 



(vr, T 



= a + hF + cF\ (3) 



resulting in a = -0.44 ± 0.06, b = 1.32 ± 0.03 GeV-^/^ and c = -0.033 ± 0.004 GeV-^. Up to the top SPS energy 
the mean pion multiplicity in p -f p interactions is approximately proportional to F . A fit of Eq. ([3]) in the range 
2 < F < 5 GeVi/2 with a = c = yielded a value of = 1.063 ± 0.003 GeY'^/^. 

For central Pb-|-Pb and Au-f Au collisions the energy dependence is more complicated as is seen in the right panel 
of Fig. [3] where the difference between the data and the p-|-p parameterization Eq. ([3]) is plotted. Below AQA GeV the 
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ratio {tt)/{Nv,) is lower in A+A collisions than in p+p interactions (pion suppression) while at higher energies this 
ratio is larger in A+A collisions than in p + p interactions (pion enhancement). The transition from pion suppression 
to pion enhancement is clearly demonstrated in the figure. A linear fit for F < 1.85 GeY^^^ using Eq. ([3]) with c = 
gave a = -0.45 ± 0.05 and b = 1.03 ± 0.05 GeV'^/'^ . The slope parameter fitted in the range F > 3.5 GeV^/'^ 
was b = 1.33 ± 0.03 where the lowest point at the top RHIC energy was excluded from the fit. Thus, in the region 
15-40A GeV between the highest AGS and the lowest SPS energy the slope increases by a factor of about 1.3. 
Fig. |4] shows the full phase space ratios (K~'')/(7r+) and £^3 = ((A) + (K + K))/(7r) as a function of collision energy 
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FIG. 4: Left: Energy dependence of the (K-'-)/{7r-'-) ratio measured in central Pb+Pb and Au+Au [isl. fl^) collisions (full 
symbols) compared to the corresponding results from p + p(p) reactions (open cirles). Right: Energy dependence of the relative 
strangeness production as measured by the Es ratio (see text) in central Pb+Pb and Au+Au collisions (full symbols) compared 
to results from p + p(p) reactions (open circles). The curves in the figures show predictions of various models described in the 
text. 

in the left and right panels, respectively. Kaons are the lightest strange hadrons and (K^) accounts for about half of 
all the anti-strange quarks produced in Pb+Pb collisions at AGS and SPS energies (a detailed explanation is given 
below) . In the Eg ratio all main carriers of strange and anti-strange quarks are included. The values for Pb+Pb and 
Au+Au collisions were calculated using data from this paper and Refs. 0, [IB, [l3j The neglected contribution 
of A and other hyperons and anti-hyperons is about 10% at SPS energies. Both the (K^)/(7r-'-) and Eg ratios are 
approximately, within 5% at SPS energies, proportional to the ratio of total multiplicity of s and s quarks to the 
multiplicity of pions. It should be noted that the {K'^)_l(Tr~^) ratio is expected to be similar (within about 10%) 
for p + p, n + p and n + n interactions at 158^1 GeV [13|, whereas the Eg ratio is independent of the isospin of 
nucleon-nucleon interactions. Calculating the Eg ratio for the p + p data Q is significantly more precise than taking 
the ratio (K+)/(7r+). 

It is seen from Fig. [4] that a steep increase of both ratios in the AGS energy region is followed by a turnover and a 
decrease around 30A GeV. The BRAHMS measurements at the top RHIG energy [Igj] indicate that the (K~'')/(7r+) 
ratio stays nearly constant starting from the top SPS energy. The RHIG results for the Es ratio are not available 
because the total multiplicity of A hyperons is not measured. For comparison the results from p + p interactions Q 
are also plotted in Fig. (4] These data show a monotonic increase with increasing energy. 

The ratios (K^) / (n^) and K^/K^ at mid-rapidity (see below for details) increase monotonically with increasing 
collision energy, as can be seen from Fig.[5l The difference between the dependence of the K"*^ and K~ yields on collision 
energy can be attributed to their different sensitivity to the baryon density. K-*- and carry a dominant fraction 
of all produced s-quarks exceeding 95% in Pb+Pb collisions at 158^4 GeV if open strangeness is considered. Because 
(K"*") = (K°) in approximately isospin symmetric collisions of heavy nuclei, the K"*" yield is nearly proportional to the 
total strangeness production and only weakly sensitive to the baryon density. As a significant fraction of s-quarks 
(about 50% in central Pb+Pb collisions at 158A GeV) is carried by hyperons, the number of produced anti-kaons 
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FIG. 5: Left: Energy dependence of the ratio (K^)/(7r^) measured in central Pb+Pb and Au+Au [IE, [3 collisions (full 
symbols). Right: Energy dependence of the ratio of mid-rapidity yields of K~ and K"*" mesons. The curves in the figures show 
model predictions described in the text. 



K and K is sensitive to both the strangeness yield and the baryon density. 

In Fig. [6] the mid-rapidity ratios K^/7r+ and K^/tt^ for central Pb+Pb and Au+Au collisions are shown as a 




FIG. 6: Energy dependence of the ratio /■k'^ (left) and K /vr (right) at mid-rapidity measured in central Pb+Pb and 
Au+Au coUisions. The NA44 data are taken from [20| . 

function of the collision energy. All NA49 mid-rapidity results presented here are determined from the fits of the 
rapidity spectra to Eq. ^ . The statistical and systematic errors were calculated taking into account all correlations 
between the fitted parameters. It is seen from Fig. [6] that the NA44 data point [i^l is consistent with the NA49 
results and that the STAR and PHENIX measurements at the top RHIC energies are in agreement with the trend 
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seen in the SPS results. Comparison with the left-hand plots of Fig. 0] and [5] shows that the energy dependence of the 
mid-rapidity kaon to pion ratios is similar to that of the corresponding ratios measured in full phase-space. 

Fig. [7] presents the energy dependence of the inverse slope parameter T of the transverse mass spectra of (left 




FIG. 7: Energy dependence of the inverse slope parameter T of the transverse mass spectra of (left) and K~ mesons (right) 
measured at mid-rapidity in central Pb-|-Pb and Au-|-Au [isl . [T6l | collisions. The K"*" slope parameters are compared to those 
from p + p(p) reactions [2j] in the left-hand plot (open circles). The curves in the right-hand plot represent predictions from 
various models described in the text. 

panel) and mesons (right panel) produced in central Pb+Pb and Au-|-Au collisions. One observes a plateau at 
SPS energies which is preceded by a steep rise of T measured at the AGS [l5| and followed by an indication of a 
further increase of the RHIC data Although the scatter of data points is large, T appears to increase smoothly 
in p -I- p(p) interactions [2l| as shown in the left panel of Fig. [71 

The transverse mass spectra of pions and protons are non-exponential such that the inverse slope parameter depends 
on the transverse mass interval used in the fit. The mean transverse mass (mx) — m provides an alternative charac- 
terization of the Tot spectra that avoids this problem, ft was calculated from multi-parameter fits to the measured 
data in the interval mx — to < 2.0 GeV. For pions either a sum of two exponentials or a power law function was used, 
which both yield a good description of the data. In case of kaons either a single exponential or a power law was used. 
The final (tot) — to values are the average of the result from both fit methods. The error is the quadratic sum of the 
statistical error and a systematic contribution that is estimated from the differences in the results obtained with the 
two different fit functions. The values of (tot) at AGS and RHIC energies were calculated from the spectra published 
in jl^, [l6|. The measurements of (tot) for p and p at the SPS were taken from jS^]. The energy dependence of (tot) 
for pions, kaons and protons is shown in Fig. [5] The results show that the approximate energy independence of (tot) 
in the SPS energy range is a common feature for all particles investigated. 

In conclusion, rapid changes in the energy dependence of pion and kaon production properties are observed which 
all seem to coincide in the low SPS energy range of 20-30^ GeV. This suggests that a common underlying physics 
process is responsible for these changes. 



DISCUSSION OF MODEL EXPLANATIONS 



In this section the energy dependence of pion and strangeness production properties will be discussed within various 
approaches to nucleus- nucleus collisions and compared with published model predictions. 

It was suggested in Q that a transition to a deconfined state of matter may cause anomalies in the energy 
dependence of pion and strangeness production in nucleus-nucleus collisions. This led to the formulation of the 
Statistical Model of the Early Stage (SMES) d, @| which is based on the assumption that the system created at the 
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FIG. 8: Energy dependence of the mean transverse mass (mx) measured at mid-rapidity in central Pb+Pb an Au+Au collisions 
for TT* (left) and (middle). For completeness, the previously published results on p and p ^22l\ are shown in the right panel. 
In the plots, positively (negatively) charged hadrons are indicated by the full (open) symbols. The results from AGS and RHIC 
are taken from p^. ITq]. 



early stage (be it confined matter or a QGP) is in equilibrium and that a transition from a reaction with purely 
confined matter to a reaction with a QGP at the early stage occurs when the transition temperature Tc is reached. 
For Tc values of 170-200 MeV the transition region extends from 15 A to 60 A GeV [6]. Assuming the generahzed 
Fermi-Landau conditions 0, [13, [13] for the early stage of nucleus-nucleus collisions and a proportionality of the pion 
multiplicity to the early stage entropy, the ratio {t:)/{N^) increases linearly with the Fermi measure F outside the 
transition region. The slope parameter is proportional to g^^'^ 5], where g is the effective number of internal degrees 
of freedom at the early stage. In the transition region a steepening of the pion energy dependence is due to the 
activation of a large number of partonic degrees of freedom. This is, in fact, observed in the data on central Pb+Pb 
and Au-|-Au collisions, where the steepening starts at about 20^1 GeV, as shown in Fig. [3l The linear dependence of 
(7r)/(iV„) on F is approximately obeyed by the data at lower and higher energies (including RHIC). An increase of 
the slope by a factor of about 1.3 is measured, which corresponds to an increase of the effective number of internal 
degrees of freedom by a factor of 1.3^ ^ 3, within the SMES [|. 

The (K~'')/(7r+) and Es ratios are roughly proportional to the total strangeness to entropy ratio which in the SMES 
model is assumed to be preserved from the early stage till freeze-out. At low collision energies the strangeness to 
entropy ratio increases steeply with collision energy, due to the low temperature at the early stage (T < Tc) and 
the high mass of the strangeness carriers in the confined state (the kaon mass, for instance, is 500 MeV). When 
the transition to a QGP is crossed (T > Tc), the mass of the strangeness carriers is significantly reduced to the 
strange quark mass of about 100 MeV. Due to the low mass m < T, the strangeness yield becomes (approximately) 
proportional to the entropy, and the strangeness to entropy (or pion) ratio is independent of energy. This leads to 
a decrease in the energy dependence from the larger value for confined matter at Tc to the QGP value. Thus the 
measured non-monotonic energy dependence of the strangeness to entropy ratio is followed by a saturation at the 
QGP value. Such anomalous energy dependence can indeed be seen in Fig. 2] and is, within the SMES, a direct 
consequence of the onset of deconfinement taking place at about 30A GeV. 

In the mixed phase region the early stage pressure and temperature are independent of the energy density [23 |. 
Consequently, within the SMES model this should lead to the weakening of the increase with energy of the inverse 
slope parameter T or, equivalently, the mean transverse mass (mx) in the SPS energy range (25l |. This qualitative 
prediction is confirmed by the results shown in Figs. [7] and [51 Moreover, recent hydrodynamic calculations [2^ 
that model both the deconfined and hadronic phases provide a quantitative description of the data as shown by the 
dashed-dotted curve in Fig. [71 

Several other analyses of the energy dependence of hadron production properties in central Pb+Pb and Au+Au 
collisions within various theoretical approaches support the hypothesis that the onset of deconfinement is located at 
the low SPS energies. In particular such a result was obtained from studies of hadron yields within a non-equilibrium 
hadron gas model [27| and using the momentum integrated Boltzmann equation for a description of the time evolution 
of the relative strangeness yield [28!| . Furthermore it was deduced from the experimentally measured rapidity spectra 
that within Landau's hydrodynamical model the sound velocity at the early stage of the reaction has a minimum at 
about 30 A GeV [29|. A minimum of the sound velocity is expected to occur in the phase transition domain. Moreover, 
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a simultaneous analysis of the two-pion correlation function and the transverse mass spectra found a plateau in the 
averaged phase-space density at SPS energies which may be associated with the onset of deconfinement [30| . 

Numerous models have been developed to explain hadron production in reactions of heavy nuclei without explicitly 
invoking a transient QGP phase. The simplest one is the statistical hadron gas model 3l| where independent of 
the collision energy the hadrochemical freeze-out creates a hadron gas in equilibrium. The temperature, the baryon 
chemical potential and the hadronization volume are free parameters of the model and are fitted to the data at each 
energy. In this formulation, the hadron gas model cannot predict the energy dependence of hadron production so 
that an extension of the model was proposed, in which the values of the temperature and baryon chemical potential 
evolve smoothly with collision energy [32] ■ The energy dependence of the Es ratio calculated within this extended 
hadron gas model is compared to the experimental results in Fig. (4] By construction, the prevailing trend in the 
data is reproduced by the model but the decrease of the ratio between 30^ and 80A GeV is not well described. The 
measured strangeness to pion yield in central Pb-|-Pb collisions at 158A GeV is about 25% lower than the expectation 
for the fully equilibrated hadron gas 32, [s^. Obviously the non-equilibrium hadron gas models 27, 34, 35 1 with 
the parameters fitted separately to the data at each energy describe the experimental results significantly better. 
Interestingly, it was found in [36| that the transition from baryon to meson dominated freeze-out conditions happens 
to be located at low SPS energies. 

Dynamical models of A-l-A collisions, such as RQMD UrQMD [3§| and HSD [i^l treat the initial nucleon- 
nucleon interactions within a string-hadronic framework. In addition these models include effects such as string-string 
interactions and hadronic re-scattering which are expected to be relevant in A-f A collisions. The predictions of the 
RQMD [33,[3i, UrQMD SO, 4a,ii] and HSD [4l| models are shown in Figs. El HEl [6] and El It is seen that aU these 
models, like the hadron gas model, fail to describe the rapid change of the hadron production properties with collision 
energy in the low SPS energy range. It was recently shown that the maximum in relative strangeness production 
can be reproduced by invoking an unusually long lifetime of the fireball at low SPS energies which decreases with 
the collision energy [43|. This assumption is however difficult to justify with the dynamical models of the collision 
process [s^, 42 [ and the results on the energy dependence of the two-pion correlation function 44 1. Finally, the 



onset of the step-like structure in the energy dependence of the inverse slope parameter of the mx spectra can be 
reproduced within the hydrodynamical model by introduction of a rapid change of the freeze-out conditions at low 
SPS energies [i^. However, this assumption does not explain the increase of the T-parameter suggested by the 
RHIC results. Thus, one can conclude that the models which do not invoke the onset of deconfinement at the low 
SPS energies can not explain the energy dependence of hadron production properties in central Pb+Pb (Au+Au) 
collisions. 



SUMMARY 



In summary, new results on charged pion and kaon production in central Pb-|-Pb collisions at 20A and 30^ GeV 
were presented and compared to measurements at lower and higher energies. A change of energy dependence is 
observed around 30A GeV for the yields of pious and kaons as well as for the shape of the transverse mass spectra. 
Available model explanations are discussed. At present a reaction scenario with the onset of deconfinement at low 
SPS energies best reproduces the data. 
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